CP asymmetry in B ^ (j)Ks in a SUSY SO (10) GUT 
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We study the B 4>Ks decay in a SUSY SO(IO) GUT. We calculate the mass spectrum of 
sparticles for a given set of parameters at the GUT scale. We complete the calculations of the 
Wilson coefficients of operators including the new operators which are induced by NHB penguins at 
LO using the MIA with double insertions. It is shown that the recent experimental results on the 
time-dependent CP asymmetry S^k in B — > 4>Ks, which is negative and can not be explained in SM, 
can be explained in the model where there are flavor non-diagonal right-handed down squark mass 
matrix elements of 2nd and 3rd generations whose size satisfies all relevant constraints from known 
experiments (t — > /^7, B — > Xsy,Bs fj,'^n~,B — > Xsfi'^ fi~ , B — > Xsg,AMs, etc.). At the same 
time, the branching ratio for the decay can also be in agreement with experimental measurements. 



I. INTRODUCTION 



Great progresses have been made on the flavor physics in recent years. Among them the progress in neutrino 
' physics is particularly impressive. Atmospheric neutrino 0] and solar neutrino U experiments together with the 
[ reactor neutrino 0, Q experiments have established the oscillation solution to the solar and atmospheric neutrino 
•- anomalies, which signals the existence of new physics beyond the SM. Experiment results indicate smallness of the 
^ ' masses of neutrinos and the bilarge mixing pattern among the three generations of neutrinos. Because the small 
: quark mixing in the CKM matrix is related to the large quark mass hierarchy understanding the bilarge mixing 
pattern is somewhat of a challenge. However, if allowing asymmetric form for the mass matrix which for example may 
well be generated by the elegant Froggatt-Nielsen (FN) mechanism |^ , we can accommodate the large mass hierarchy 
' with large mixing Therefore, if one works with an effective theory, e.g., the minimal supersymmetric Standard 
I Model together with right-hand neutrinos (MSSM-I-N), at a low energy scale (say, the electro- weak scale), one can 
> content oneself by using the WWZF scenario to understand the smallness of quark mixing, and the see-saw and FN 
mechanism to understand the smallness of the masses of neutrinos and the largeness of neutrino mixing. However, if 
one works with a theory, e.g., a grand unification theory (GUT), in which quarks and leptons are in a GUT multiplet, 
one has to answer: can we explain simultaneously the smallness of quark mixing and the largeness of neutrino mixing 
Q-i in the theory? If we can, then what are the phenomenological consequences in the theory? There are several recent 
^ works to tackle these problems in SU(5), flipped SU(5), or SO(IO) GUTs [lllli|ll|ll|ll|ll|l3, which brings 
the study of GUT to a more realistic level. 

The smallest grand unification group that incorporates righthanded neutrinos required for the seesaw mechanism, 
is SO(IO). It has been shown that the observed bilarge mixing naturally leads to large flavor non-diagonal down-type 
squark mass matrix elements of 2nd and 3rd generations in a SUSY SO(IO) llj. 
' The measurements of the time dependent CP asymmetry Sj/^k in _B ^ J/ipKs have established the presentence 
of CP violation in neutral B meson decays and the measured value [T6| 

Sj/^K = sin(2/3( J/V;i^s))worid-avc = 0.734 ± 0.054. (1) 

is in agreement with the prediction in the standard model (SM). Recently, various measurements of CP violation in 
B factory experiments have attracted much interest. Among them [l7lll8| . 



S^Ks = -0.39 ± 0.41, 2002 World - average 

S<t,Ks = -0.15 ± 0.33, 2003 World - average (2) 

is especially interesting since it deviates greatly from the SM expectation 

S^Ks = sin{2P{(bKs)) = sm{2(3{J/^Ks)) + 0{\^ (3) 

where A ~ 0.2 appears in Wolfenstein's parameterization of the CKM matrix. Though the impact of these experimental 
results on the validity of CKM and SM is currently limited by experimental uncertainties, they have attracted much 
interest in searching for new physics [millllllil and it has been shown that the deviation can be understood 
without contradicting the smallness of the SUSY effect on B — > J/tpKs in the minimal supersymmetric standard 
model (MSSM) ||23, |2j|. Motivated by SUSY GUTs, the time dependent CP asymmetry in B ^ <j)Ks has been 
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studied in the SUSY models with a large mixing between the 2nd and 3rd generation in the down-type squark mass 
matrix at the Mgugy ^ TeV scale ^25]. The asymmetry has also been examined in the SUSY SU(5) framework [26| . 
It is shown in ref. 26] that the possibility of large deviations from SM in S'^if is excluded in the case of only having 
the flavor non-diagonal down-type squark mass matrix element of 2nd and 3rd generations in the RR sector (i.e., ^23*^ 
non-zero with (5|J^^ = (A'/|^^)23/m?, where (M|^^)y is the flavor non-diagonal squared right-handed down squark 

mass matrix element and rriq is the average right-handed down-type squark mass) due to the bound on 52§^ from 
Br{T fij). In this paper we investigate the decay B (pKs in the SUSY SO(IO) framework. Our results show 
that the time dependent CP asymmetry in i? — > (/)ifs can sizably deviate from SM after imposing the constraint from 
Br{T — !■ /i7), in contrast with the claim in the literature, because we have included the double insertion contributions 
in penguin diagrams for the relevant Wilson coefflcients. Furthermore, the contributions from neutral Higgs penguins, 
which are not considered in ref. 26], have been included in this paper. 

We need to have new CP violation sources in addition to that of CKM matrix in order to explain the deviations 
of Sc/iKs from SM. There are new sources of flavor and CP violation in SUSY SO(IO) GUTs. In such kind of models 
there is a complex flavor non-diagonal down-type squark mass matrix element of 2nd and 3rd generations of order 
one at the GUT scale (ll| which can induce large flavor off-diagonal couplings such as the coupling of gluino to the 
quark and squark which belong to different generations. These couplings are in general complex and consequently can 
induce CP violation in flavor changing neutral currents (FCNC). It is well-known that the effects of the counterparts 
of usual chromo-magnetic and electro-magnetic dipole moment operators with opposite chirality are suppressed by 
ms/mi, and consequently negligible in SM. However, in SUSY SO(IO) GUTs their effects can be significant, since 
62^^ can be as large as 0.5 11]. 

For the b s transition, besides the SM contribution, there are mainly two new contributions arising from the 
QCD and chromo-magnetic penguins and neutral Higgs boson (NHB) penguins with the sparticles propagating in the 
loop in SUSY models. The contributions to the relevant Wilson coeflScients at the mw scale from the diagram with 
the gluino propagating in the loop have been calculated by using the mass insertion approximation (MIA) l3C| with 
double insertions in ref. p3 ] . We calculate the contributions from the diagram with the chargino propagating in the 
loop using MIA with double insertions because chargino contributions can be signiflcant when the left-right mixing 
between squark masses is large. 

As it is shown that both Br and CP asymmetries depend significantly on how to calculate hadronic matrix elements 
of local operators 24]. Recently, two groups, Li et al. [s^ [s^ and BBNS [33. IssI ]. have made significant progress 
in calculating hadronic matrix elements of local operators relevant to charmless two-body nonleptonic decays of B 
mesons in the PQCD framework. The key point to apply PQCD is to prove that the factorization, the separation 
of the short-distance dynamics and long-distance dynamics, can be performed for those hadronic matrix elements. 
It has been shown that in the heavy quark limit (i.e., mi, — > 00) such a separation is indeed valid and hadronic 
matrix elements can be expanded in as such that the tree level (i.e., the a° order) is the same as that in the naive 
factorization and the as corrections can be systematically calculated |34j. Comparing with the naive factorization, 
to include the as correction decreases significantly the hadronic uncertainties. In particular, the matrix elements of 
the chromomagnetic-dipole operators Qg^ have large uncertainties in the naive factorization calculation which lead to 
the significant uncertainty of the time dependent CP asymmetry in SUSY models [2^. The uncertainties are greatly 
decreased in BBNS approach [ssj. The hadronic matrix element of operators relevant to the decays B 4>Ks up to 
the as order have been calculated in BBNS approach in ref. ^ . 

Using the BBNS approach to calculate hadronic elements to the as order and Wilson coefficients in MIA with 
double insertions, we show in this paper that in the SUSY 5*0(10) framework in the reasonable region of parameters 
where the constraints from r ^7, Bs — Bs mixing , r(6 57), T{b ~> sg), T{b — > sii^ii~), and B are 
satisfied, the branching ratio of the decay for B — > 4)Ks can be smaller than 1.6 x 10~^, and the theoretical prediction 
for Sif)K can be in agreement with the data in la experimental bounds and even can be as low as —0.6. 

The paper is organized as follows. In Section II we describe the SO(IO) models we used in the paper briefly. In 
section III we calculate sparticle spectrum using revised ISAJET. In section IV we give the effective Hamiltonian 
responsible for B 4>Ks in the model. In particular, we give the Wilson coefficients of operators using MIA with 
double insertions. The Section V is devoted to numerical results of the time dependent CP asymmetry and branching 
ratio for the decay B —> 4>Ks. We draw conclusions and discussions in Section VI. 

II. NEUTRINO BILARGE MIXING AND LARGE B-S TRANSITIONS IN SO(10)GUT 

It is well-known for a long time that SO(IO) GUTs can naturally incorporate the seesaw mechanism (s^ which 
is the simplest way to understand small neutrino masses. Recently a number of SUSY SO(IO) models which can 
explain neutrino data and quark mixing have been proposed 0, llll fl^ and some phenomenological consequences 
of the models have been analyzed 0, 0, O] . It was first pointed out in that the induced large mixing in the 
down-type squark mass matrix in SUSY 5*0(10) models have interesting consequences in low energy B physics. For 
specific, we use the model in ref. IT] and review the main points of the model. The details of the model can be found 
in Ref. [l^. 
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In order to accommodate both CKM mixing among quarks and MNS mixing among leptons, one needs to have an 
asymmetric down- type Yukawa matrix Yf^. So instead of the usual superpotential |3^ . one assumes 

W = i(r„)yl6a6jl0, + hYd)^Jl6^16,^Wd■ (4) 
Z Z -IVlpi 

Because of the combination of the Higgs multiplet 45, whose VEV (45) 7^ breaks SO(IO), and the Higgs in 10, the 
effective Yukawa couphng being either in 10 (symmetric between two 16's) or 120 (anti-symmetric between two 16's) 
representations, the matrix Yii can now have a mixed symmetry. Setting the breaking chain to be S'O(IO) SU{5) — > 
SM, we have the Yukawa couphngs in the MSSM-I-N (the MSSM with right-handed neutrinos) as 

W = {Y^\Q;U,H^ + {Y^%UN,Hu 

+ {V*Y^U),,Q,D^Ha + {V*Y^U\,E,LjHa + ]^M,,N,Nj, (5) 

where Yf (i=u, d) is the diagonal and we have absorbed the phase matrix Ql into Qi and E'i, the phase matrix Qn into 
{UD) multiplet or {EV*) multiplet, and the phase matrix 8^ into the Majorana mass matrix M, respectively. Here V 
and U are the usual CKM mixing matrix and neutrino mixing matrix (MNS matrix) respectively. We have taken Yu 
and M simultaneously diagonal. Such a situation could result from simple U{1) family symmetries and in a 5*0(10) 
with hierarchical Kf* and right-handed neutrino masses, the choice of having such simultaneous diagonalization looks 
rather plausible |39j . 

After integrating out the right-handed neutrinos in Eq. Q, the light neutrino masses are determined from the 
superpotential 

W = \{Y^UM-%{Y^),{UH,,){L,H^), (6) 

which leads to the light Majorana neutrino mass matrix {m^)nm = {Y^^ /'2AIi)U*j^e~'^^^U*^j^, where e^^^ is the phase 
of the diagonal element Mi of the diagonal matrix M. Therefore, the two mass splitting data and bilargc mixing of 
neutrinos can be explained. 

The Yukawa coupling of "third-generation" neutrino is unified with the large top Yukawa coupling due to the 
50(10) unification. Nevertheless, this "third-generation" neutrino is actually a near-maximal mixture of t'^j and Vt, 
which comes from the large mixing angle in atmospheric neutrino oscillation. Because the "third-generation" charged 
leptons, down-type quarks and neutrinos are in a SU{5) multiplet, the SU{5) multiplet with the large top Yukawa 
coupling contains approximately 

5; = 5; cos 5* sine, (7) 
where 9 ~ 45° is the atmospheric neutrino mixing angle, and 

5; = (6^6^5^^.„r) (8) 

5; = (s^s^s^^.^,^i). (9) 

The similar mixing exists in right-handed b and s quarks. However, we do not worry about its effects because there 
is no charged-current weak interaction on right-handed quarks. Therefore the mixing among right-handed quarks 
decouples from low-energy physics. 

We now come to the main outcome in the framework which has phenomenological consequences in low-energy B 
physics. Because one works in SUSY models there is the corresponding mixing among squarks which would yield 
observable effects. The top Yukawa coupling generates an 0(1) radiative correction to the mass of ssin0 -|- bcosd, 
which leads to a large mixing between s and b at low energies. This large mixing in turn generates interesting effects 
in i?-physics. 

To be specific, assuming that at the scale which is above the GUT unification scale, say, near the Planck scale, 
one has the universal soft terms: scalar mass mo, gaugino mass M1/2, trilinear and bilinear couplings Aq and Bq, 
then large neutrino Yukawa couplings involved in the neutrino Dirac masses can induce large off-diagonal mixing in 
the right-handed down squark mass matrix through renormalization group evolution between and Mqut and the 
induced mixings will generally be complex with new CP violating phases. The induced off-diagonal elements in the 
mass matrix of the right-handed down squarks d^ are given by (in the basis in which Yjj is diagonal) 



^ -T J nm 



where M5 is the SU{5) breaking scale and 
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FIG. 1: The mass spectrum versus mo for fixed Mi/2=200GeV, tan/3=40, Si^nn = (0.12 + 0.03i), and sign(^)=+l without 
the constraints from the low energy experiments imposed, (a) is for Ao — 0. (b) is for ylo=-1000GeV. 



where e"^"" is the phase from (8fl)„„. Note that these phases are not relevant to any other low energy physics. 
Substituting C/23,33 into ea. l(TT|l . one obtains 

[y«ty«]23 = 0.5e-'('^^^'^'^3^')(mt(MG)/178 GcV)^ (12) 

where rntiMc) is the top quark mass at Mq- From eqs. H1UI12|I . the mass insertion parameter 62^^ can easily be of 
order one. 



III. MASS SPECTRA AND THE PERMITTED PARAMETER SPACE 



To see the impact of the induced off-diagonal elements in the mass matrix of the right-handed down-type squarks on 
B rare decays and simplify the analysis, we assume that at the GUT scale (Mg) all sfermion mass matrices except the 
right-handed down-type squark mass matrix are flavor diagonal and all diagonal elements are approximately universal 
and equal to TOq. The 2-3 matrix element of right-handed down- type squark mass matrix is parameterized by ^23*''* = 
(M|^^)23/mQ which can be treated as a free parameter of order one, as discussed in last section. Furthermore, we have 
a universal gaugino mass M1/2, a universal trilinear coupling Aq and a universal bilinear coupling Bq at Mq. With the 
requirement of radiative electro-weak (EW) symmetry breaking, we have five parameters (mp, Af]^/2i ^Oi (^j-f^, ton/?) 
plus the sign of fi as the initial conditions for solving the renormalization group equations (RGBs). 

We require that the lightest neutralino be the lightest supersymmetric particle (LSP) and use several experimental 
limits to constraint the parameter space, including l)the width of the decay Z XiXi is less than 4.3 MeV, and 
branching ratios of Z ^ X1X2 ^ ^ X2X2 ^-^^ less than 1 x 10~^, where Xi is the lightest neutralino and X2 is the 
other neutralino, 2) the mass of light neutral even Higgs can not be lower than 111 GeV as the present experments 
required, 3) the mass of lighter chargino must be larger than 94 GeV as given by the Particle Data Group [s^, 4) 
sneutrinos are larger than 94GeV, 5) seletrons are larger than 73GeV, 6) smuons larger than 94 GeV, 7) staus larger 
than 81.9 GeV. 

We use revised ISAJET to do numerical calculations. We find that the parameter (-^|^^)23 does not receive any 
significant correction and the diagonal entries are significantly corrected, which is in agreement with the results in 
ref. [25. We scan mo, M1/2 in the range (100, 800) GeV for given values of Ao,tan/3 and sign(^)=+l* We impose 
the constraints from the relevant low energy experiments such as i? — > Xsj, etc. (for the detailed analysis of the 
constraints, see section V). 

As an illustration, the mass spectra without and with the constraints from the low energy experiments are given 
in Fig. 1 and 2, respectively, where (a) and (b) are for Aq = 0, — lOOOGeV respectively. Fig. 1 is plotted for 
Ml/2 = 200G'ey,tan^ = 40 and Fig. 2 for M1/2 = SOOGeF, tan/3 = 40. One can see from the Figs. [2 H that 
the mass spectrum lifts when |^o| decreases. When the constraints from the low energy experiments are imposed, 
some mass spectra are excluded and for the allowed spectra the masses of sparticles are similar to those without the 
constraints. 



* In the case of sign(/^)= -1, the constraint from B X^f on the parameter space is too stringent, in particular, for large tan/3 I56tl57l . 
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FIG. 2: The mass spectrum versus mo for fixed Mi/2=500GeV, tan/3=40, 523flfl = (0-12 + 0.03i), and sign(^)=+l with the 
constraints from the low energy experiments imposed, (a) is for Aq — 0. (b) is for ylo=-1000GeV. 



IV. EFFECTIVE HAMILTONIAN FOR s TRANSITION 



The effective Hamiltonian for 6 ^ s transition can be expressed as[2^ |4| 



+Cr-, Q77 + Csg Qsg + C'j^ + C'^g Q'sg j + h.c. (13) 
Here Qi are quark and gluon operators and are given by ^ 

Ql = {SaPl3)v-A{Pf}ba)v-A, Q2 = {SaPa)v -A{ppbi3)v -A, 

Q3(5) = {Saba)v-A^^iqpqi:i)v-(+)A, Q4(6) = iSabi3)v-A^^{qpqa)v-{+)A, 

Q q 

Q7{9) = ■:^^Saba)v-A'^eq{qi3qi3)v + (-)A, ^8(10) = ^ (ga^/j)y-A eg(g;3ga) 

? 9 



(3ll(l3) = {sb)s+py^ —{qq)s-i+)p, 
g 

Ql2(14) = (Si^j)s+PE qi)s-{+)P : 

Till) 

q 

Qi5 =sa^^(l + 75)6V ^ga^.(l+75)g, 



Qie = Si(T^''(l + 75) bj — gj- (T^^(l + 75) , 

(377 = ^TO&SqO''^''-F}'I.(1 +75)^/3, 

= |jjmbS„a'^''G;,^(l + ^,)bp, (14) 

where (gig2)v±yi = qiY{^ ± 75)92, {qiq2)s±p = ± 75)'?2 p = u,c, g = u,d,s,c,b, Cq is the electric charge 
number of q quark, Aq is the color SU(3) GcU-Mann matrix, a and (3 are color indices, and i^^jy {Gp^u) are the photon 



t For the operators in SM we use the conventions in Ref. l35l where Qi and Q2 are exchanged each other with respect to the convention 
in most of papers. 

t Strictly speaking, the sum over q in expressions of Qi (i=ll,...,16) should be separated into two parts: one is for q=u, c, i.e., upper type 
quarks, the other for q=d, s, b, i.e., down type quarks, because the couplings of upper type quarks to NHBs are different from those 
of down type quarks. In the case of large tan/3 the former is suppressed by tan~^ (3 with respect to the latter and consequently can be 
neglected. Hereafter we use, e.g., CJ^j to denote the Wilson coefficient of the operator Q\\ = (sb)g^p |^(cc)s_p. 
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(gluon) fields strength. 

The primed operators, the counterpart of the unprimed operators, are obtained by replacing the chiralities in the 
corresponding unprimed operators with opposite ones. The SUSY contributions to Wilson coefRcients have been 
calculated by using the vertex method in ref . |29j | . The SUSY contributions due to gluino box and penguin diagrams 
to the relevant Wilson coefficients at the mw scale in MIA with double insertions, as investigated in ref which are 
non-negligible if the mixing between left-handed and right-handed sbottoms is large have been given in refs.pfll l4fil|. 
We calculate the chargino contributions in MIA with double insertions and results are 



Qas 
1,2 



with C^'^ J 

Co, — 



~2V^lV:^htFe{x)5l!^''] 

^32 

-2VaV,lhtFeix)S!^^ 
Cs = C3 

Cq = C4 

Cr, - ^^[^T^iF;i(^^i2(x)<52"3"+F{2(x)<5^f«<5^f«) 

-^VaVMF^,ix)s^i:^ + F[,ix)siii:H-i:^) 

^32 

-2V^lV*2htFl2{x)^^i'^] 

+ l-^[-^ynU:,K{F,,{x)6-i^ + ^^^4(^)^2t''^3t'') 
9 mbm.± A32 

uLRi 



uLR-i 



+ V,2U:2hthf,F3i{x)d^i 

Cs, = ^^[^VaV:,{F2{x)6^,^^+F^ix)6^ii^^S};i:^) 
24 m^± A32 

~^VaVMF2{x)Sl^i''^ + F!,{x)5^^H^^^) 

^32 

-2Vi^V:2htF2{x)5lt''] 

1 K* 

' ''V,iU*,h,{F,{x)8-^^+F',{x)5-,^^5-,^^) 



6 mbm^± K^2 



uLR} 



'V,2U*2hthbF^{x)5^^ 

167r^ mi 

=2 

IGtt^ mi 

= 0, i = 12,14,15,16 (15) 



^11 - T7^^ — i^Qi T^qJ 

13 1 R_2 m, ^1 V2J' 



1 mim-± ^2 I 2 T^* 

TTY 2 2 [--f7rV^iU^2'T'b[^bO[X)d23 + ^bo(^)(>23 '^33 ) 

V^2U*2hthbFbo{x)5^^'^'] 



5 The operator ofQ^'l, is defined as Qi = ^ [5(1 + 75 [/I], Q'^ = ^ [5(1 -75 [/I], = ^[s{l + 'y5)b][h5l], Q'^ = ^ -75)b] f75i] 
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Cn^ — 



I mim~± 



2sl 



'■HO 



^32 



LL 
23 



uLR-\ 



(16) 



where ht = —m^ — , ^ifc = nT"^ — , '"s = — ^» = — r^, and x = m?/m?± with ma and m-± being the common 

squark mass and chargino masses respectively. The repeating indices i should sum over from 1 to 2. The one-loop 
functions in Eq. (|15|) and p6|) are given i n Ap pendix A. We have checked that the Wilson coefficient Cr^ with single 
insertion is the same as that given in Ref. 42j. Differed from the single insertion results, the LR or RL insertion also 
generates the QCD penguin operators when one includes the double insertions. 

For the processes we are interested in this paper, the Wilson coefficients should run to the scale of 0{mi,). Ci — Cio 
are expanded to 0{as) and NLO renormalization group equations (RGEs) should be used. However for the Csg and 
Qlj, LO results should be sufficient. The details of the running of these Wilson coefficients can be found in Ref. 
pO| . The one loop anomalous dimension matrices of the NHB induced operators can be divided into two distangled 
groups ^3 



and 
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Qii Qi2 


Qii 


-16 


Ql2 


-6 2 



(17) 





Ql3 


Qii 


Ql5 


Ql6 


Ql3 


-16 





1/3 


-1 


Ql4 


-6 


2 


-1/2 


-7/6 


Ql5 


16 


-48 


16/3 





Ql6 


-24 


-56 


6 


-38/3 



(18) 



Here and hereafter the factor, is suppressed (i.e., the anomalous dimension matrix for (5ii.i2 is fj etc.) 
For Q'^ operators we have 



and 
(') ;-i 



(19) 



Because at present no NLO Wilson coefficients , 1=11,.. .,16, are available we use the LO running of them in the 
paper. 

There is the mixing of the new operators induced by NHBs with the operators in SM. The leading order anomalous 
dimensions have been given in Rcfs.[3.l45j|. We list those relevant to our calculations in the following. Defining 



O, = 



9 



167r2 



one has 



(RD) _ 
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The mixing of Qii,i2 onto the QCD penguin operators is 



12+j, i = 


1,2, 








Oi 


-1/3 


1 


O2 


-1 





O3 


28/3 


-4 


O4 


20/3 


-8 



(20) 



(21) 



For the mixing among the primmed operators, we have 





Q3 


Qa Q5 


Qe 


On 


1/9 


-1/3 1/9 


-1/3 


O12 












7 



{LD' 



7 



(RD) 



and 7 



(A/'Q') 



7 



(MQ) 



(22) 



(23) 



The mixing of the new operators induced by NHBs with the operators in SM has non-negligible effects on the Wilson 
coefficients of the SM operators at the 0(mf,) scale. In particular, the Wilson coefficient of the chromo-magnetic 
dipole operator C'sg at the 0{mh) scale, which has a large effect to Smk [M — 0, 77'), can significantly enhance due to 
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the mixing. To see it explicitly we concentrate on the mixing of Oi (for its definition, see Ea. 120(1 ') onto Qgg. Solving 
RGEs, we have 



c=l,...A 

^(Mo) = E lalVac^VcbCb{li.o)l{lcc- l&g&g) 



a,b=l,...,4 

where V and 7aa are given by 

+ 2/3o /) V-^ = diag(7ii,722,733,744). 

with / being the 4x4 unit matrix. Using 

Cairn) = Ci{llo)5al 

and Eq. (|2U, Eq. (|23I) reduces to 

CM = 0.68C8<,(/io) - 3.2Ci3(Aio), 



(24) 

(25) 
(26) 

(27) 

(28) 

(29) 



where Ci{fio) = ^Cisifio) has been used. 

In our numerical calculations we neglect the contributions of EW penguin operators Q7,...io since they are small 
compared with those of other operators. 

The hadronic matrix elements of operators have been calculated in BBNS approach in refs. [23.135114^ . We use the 
results in ref.^^l- The effective Hamiltonian ((13(1 results the following decay amplitude for (t>Ks- 



A{B ^ cj^Ks) - ^A, 
A = A" + A°'+A", 



(30) 



where we have divided the decay amplitude into three parts in order to see explicitly the effects of new operators in 
the SUSY SO(IO). One, A", has the same form as that in SM, the second, A° is for primmed counterparts of the 
SM operators, and the third, A", is new which comes from the contributions of Higgs penguin induced operators. In 



A" = {cly\sj^s\0){K\srb\B)x ^ VptV;, 

p—u,c 

A°' = A°{C, ^ C[), 



03 + + as - 2 («7 + 09 + afg) 



(31) 



(32) 



where Oi's have been given in Refs.|2Ct l47j. The hadronic matrix element of the vector current can be parameterized 
as (if|s7''&|B) = Ff^^{q^){p% +p^) + {Ff^^^{q^) - Ff^^{q^)){m% - m]i)qi'/q^. For the matrix element of the 
vector current between the vacuum and (/>, we have ((/)|s7p6|0) — ms fibf-t - 
A" in eg. 1(30(1 . to the as order, in the heavy quark limit is given as|24| 



A" = A"(C,) + A"(C, ^CO, 
A"(C,) = {<t>\-s^^s\Q){K\-srh\B) {-VaV:, 



04 + 



nib 



Arris 



:ai2 



vrib 



■ ai5 



(33) 



where a^'s have been given in Refs. |24l |46| . 

The time-dependent CP-asymmetry S^k is given by 



2ImA 



4>K 



Here is defined as 



q\ A{B^4>Ks) 



pJsAiB-.cj)Ks)' 



(34) 



(35) 
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V. NUMERICAL RESULTS 
A. Parameters input 

In our numerical calculations the following values are needed: 

• Lifetime, mass and decay constants 

t{B°) = 1.56 X 10^i2g, Mb = 5.28GeV, = 4.2GeV, 

TOc = l.SGeV, rus = lOOMeV, = 0.237GeV, 

• Wolfenstein parameters 

We use the Wolfenstein parameters fitted by Ciuchini et al|5^: 

^ = 0.819 ±0.040 A = 0.2237 ±0.0033, 

p = p(l - aV2) = 0.224 ± 0.038, p = 0.230 ± 0.039, 

fj = r]{l- = 0.317 ± 0.040, = 0.325 ± 0.039, 

7 = (54.8 ±6.2)°, y^p^+v^ = 0.398 ±0.040. 



Form factors 

In the paper we need the form factors: F^^^'^ (0) = 0.34. 



(36) 



(37) 



B. Constraints from experiments 

We impose two important constraints from B X^j and Eg —^ fj,~^ fi~ . Considering the theoretical uncertainties, we 
take 2.0 x 10~^ < Bi{B Xg'y) < 4.5 x 10~^, as generally analyzed in literatures. Phcnomcnologically, Br{B Xg'y) 
directly constrains \Crj{mb)\'^ + \CY^{mi,)\'^ at the leading order. Due to the strong enhancement factor mg/mi, 

associated with single 62^^^^^^ insertion term in Cj^{mi,), 62^^^^^^ {''^ 10~^) are more severely constrained than 
^23'^^^^^- However, if the left-right mixing of scalar bottom quark S^^^'"^^^ is large (~ 0.5), ^23'^^^^^ is constrained 
to be order of 10~^ since the double insertion term jg ^^jgQ enhanced by rrig/mb. Nevertheless, 

in the large tan/3 case the chargino contribution can destructively interfere with the SM (plus the charged Higgs) 
contribution so that the constraint can be easily satisfied. The branching ratio Bs — > p.^ p,^ in SUSY models is given 
as 



BT{B,~.p+p-) - -g^ml^TsjIJA.lVl - 4m2[(l - Arf^)\CQMb) ~ C'q^M'' + 

|CQ,(mf,) - C^,(TOf,) + 2m(C7io(mb) ~ C[^{mkW] (38) 

where fh — nii^JmB^- In the middle and large tan/3 case the term proportional to (C'lg — CJq) in Eq. H38() can be 
neglected. The new CDF experimental upper bound of Br(i?s p^ p~) is 5.8 x 10~^ |^ at 90% confidence level. 

To translate it into the constraint on , we have 

yii,i3 ' 

^J\CQ,A^w) ~ C'QjmwW + \Cq,A^w) - C'^JmwW < 0-047 (39) 

Because the bound constrains |Cq. — Cg. | (i=l, 2), ^ and there is a flavor non-diagonal element in the right-hand 
down-type squark mass matrix at the high scale in the SUSY SO(IO) model, we could have values of |Cq. | and \C'q.\ 
larger than those in constrained MSSM (CMSSM) with universal boundary conditions at the high scale and scenarios 
of the extended minimal flavor violation in MSSM (2^] in which |Cq. | is much smaller than |CqJ. However, since the 
flavor non-diagonal element in the right-hand down-type squark mass matrix at the high scale has almost no impact 
on the decoupled limit of Higgs boson masses (i.e., to^o ~ to^o ~ rnjf±, sin(Q; — /3) ~ 1), the decoupled Higgs sector 
in CMSSM remains in our model unless the non universal conditions at the high scale are assumed j^j. Therefore, 



^ Cq' ^ are the Wilson coefficients of the operators Qi \ which are Higgs penguin induced in leptonic and semileptonic B decays and their 
definition can be found in Ref. ISTl . By substituting the quark-Higgs vertex for the lepton-Higgs vertex it is straightforward to obtain 
Wilson coefficients relevant to hadronic B decays. 
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FIG. 3: The correlation between S^Kg and Br(_B <i>Ks) for Aq — 0. (a) is for the SUSY contributions with only gluino 
propagated in the loop, (b) is for the all contributions included. Current la bounds are shown by the dashed lines. 



Cq^ ~ "F^q! ^-"^ '^^^ model, which leads to Cq|^ being almost zero, and consequently eq. (|39|l give a stringent 

constraint on Cg^^. At the same time we require that predicted Br of i? ^ X^/^^/i^ falls within 1 a experimental 
bounds. 

We also impose the current experimental lower bound AMs > 14.4ps~^ p3 | and experimental upper bound Br(i? — s- 
Xsg) < 9% [s^. Because 62^^^^^^ is constrained to be order of 10"^ by Br(_B Xgj), their contribution to AMs 
is small. The dominant contribution to AAf^ comes from 82^^^^'^^ insertion with both constructive and destructive 
effects compared with the SM contribution, where the too large destructive effect is ruled out, because SM prediction 
is only slightly above the present experiment lower bound. 

As pointed out in section III, due to the gluino-sbottom loop diagram contribution and the mixing of NHB induced 
operators onto the chromomagnetic dipole operator, the Wilson coefficients Cg^^ can be large, which might lead to a 
too large Br of B — > Xsg. So we need to impose the constraint from experimental upper bound Br(i3 Xsq) < 9%. 
A numerical analysis for Cg^— has been performed in Ref . [isj . We carry out a similar analysis by setting both Csg 
and Cgg non-zero. 

Furthermore, as analyzed in ref. |26| . there is the correlation between flavor changing squark and slepton mass 
insertions in SUSY GUTs. The correlation leads to a bound on ^23^^ from the rare decay t — > ^j. We update the 
analyses with latest BELLE upper bound of Br(r /i7)< 3.1 x 10~^ at 90% confidence level. 



C. Numerical results 



In numerical analysis we fix tan/3 = 40, sign(/i)=+l, and Aq — 0, — lOOOGeV. We scan mo,Afi/2 in the range 
from lOOGeV to SOOGeV by running ISAJET. Using the sparticle mass spectrum and mixings at the EW scale, 
we calculate Wilson coefficients and consequently Br and the time dependent CP asymmetry in _B ^ fp^^S under 
additional constraints of (Jfl^^ from Br(r —> /17), which phase varies from to 27r. 

Numerical results of the correlation between S^Ks ^nd Br(_B 4>Ks) for Aq — 0, — lOOOGeV are shown in Figs. 13 
and 01 respectively, where (a) is for the SUSY contributions with only gluino propagated in the loop, (b) is for the 
all contributions (i.e., W^^,H^, chargino, gluino, neutrilino propagated in the loop** ) included. Current Icr bounds 
are shown by the dashed lines. The case (a) describes the direct consequence of existing a ^23^^. Before we discuss 
the numerical results in detail, a remark is in place. It is shown that in MSSM NHB contributions can be significant 
with all relevant experimental constraints imposed p^ . However, NHB contributions are limited to be very small in 
our case. In case of only SM and NHB contributions included, both S^pKs ^^^d Br(i3 <pKs) almost do not change, 
compared with the SM, because the Higgs sector is nearly decoupled with ttiho ~ to^o ~ mH±, sin(Q; ~ f3) ~ 1, as 
pointed out above, in the regions of the parameter space which we take in the SO(IO) model (or other constrained 

MSSM), which leads to that c['^ are almost equal to and consequently c[i are very small due to the constraint 
from Bs — > . 



* We neglect the neutrilino contributions in numerical calculations because they are small compared with other contributions. 
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FIG. 4: The correlation between S,i,Ks a-nd Br(B (j)Ks) for Ao — — lOOOGeV. The black (red) dots are allowed by all relevant 
experimental bounds including (without including) Br(fc sg). (a) is for the SUSY contributions with only gluino propagated 
in the loop, (b) is for the all contributions included. Current la bounds are shown by the dashed lines. 

We find from Fig. 0| that in both the cases (a) and (b) there are regions of parameters where Sii,Ks falls in la 
experimental bounds and Br is smaller than 1.6 x 10^^. In the case of all contributions included the region is larger 
than that in the case with the contributions from only gluino-squark in the loop included if the constraint from b ^ sg 
is not imposed. There exist two regions of the parameter space as given in Fig.^(b) for the all contributions. The 
dense black region corresponds to the heavy SUSY spectrum region with niQ, M1/2 as large as 600 ~ SOOGeV, while 
the scattered belt region denotes Afi/2 as small as 100 ^ 200GeV and toq in the region 500 ~ SOOGeV where sleptons 
are as heavy as 500 ^ SOOGeV and Br(T /17) constrains 823^^ to be order of 0.1. That Br(6 — > sg) serves as a 
strong constraint on the parameter space is shown in Fig. 01 in particular, the Fig. 0] (b) with all the contributions 
included. In the scattered belt region of the Fig. 0] (b), with M1/2 as small as 100 ^ 200GcV, the destructive 
chargino contribution to drives Cj^ smaller than Cj^, which calls for larger C^-^ of gluino contributions to 
enhance Br(6 — > 57) cx jCy^p + |C7..yp above the experimental lower bound. However large C^^ results in large 
Cgg, which raises Br(& sg) beyond the experimental upper bound. Therefore, the most of the region without 
the b sg constraint is excluded when the b sg constraint is imposed. In the case of Aq = 0, Fig. |31 (b) 
corresponds to only the heavy SUSY spectrum region because the light mass spectrum region, which corresponds to 
lOOGeT^ < mo, M1/2 > 4006*6^ roughly, is excluded by Br(i3 Xs'j) and the experimental lower bound of m/jo and 
has no results with S'^Xs near 0, while Fig.|2|(a), which corresponds to the region of the parameter space with as 
small as 100 ~ 200GeV and mo in the region 500 ^ SOOGeV, has some results with S^Ks iiear 0, because Br(6 37) 
constraints can be satisfied easily in the latter case. Comparing Fig. ^ with Fig. O the conclusion is that the SUSY 
effects on S'^/c^ for a negative large Aq (say, |^o| ^ "^0,-^^1/2) are larger than those for Aq = 0. The reason is as 
follows. One needs to have a large fj, in order to have a large RL mixing between down squarks, i.e., a large S^^^, and 
consequently a large induced <5f;f^. Due to the RGE running, a negative large ^0 can drive fi at the low energy large, 
whereas ^0 = or a small ^0 fails. And at the same time, a large /i can make the chargino contributions significant, 
which make the real part of Csg + Cg^ smaller than that in the small Aq case^^ . 

We vary tan/? and find similar results for tan/3 = 10,30. The numerical results are obtained for lOOGeT^ < 
Too, Afi/2 < SOOGeT^. For fixed nig, the Wilson coefficient Gg^ is not sensitive to the variation of the mass of squark 
in the range about from lOOGeV to 1500GeV. Therefore, the numerical results are not sensitive to mo for fixed M1/2 
and would have a sizable change when M1/2 decreases. 

VI. CONCLUSIONS AND DISCUSSIONS 

In summary we have calculated the mass spectrum and mixing of sparticles in the SUSY SO(IO) GUT. We have 
calculated the chargino contributions to Wilson coefficients at LO using the MIA with double insertions in SUSY 
models. Using the Wilson coefficients and hadronic matrix elements previously obtained, we have calculated the 



tt In the paper we do not consider the cancellation mechanism "63*1 in the analysis of electric dipole moments (EDMs) of electron and 
neutron in SUSY models and assume fi, Ao are real. 



12 



time-dependent CP asymmetries Scj,k and branching ratios for the decay B 4>Ks. It is shown that in the reasonable 
region of parameters where the constraints from r — > //7, Bg — Bg mixing , T[b — *■ 57), r(6 — > sg), T{b — > sfi'^fi~), 
and B /^^M^ are satisfied, the branching ratio of the decay for B <t>Ks can be smaller than 1.6 x 10~^, and 
S^iKs can be negative. In some regions of parameters S^Ks can be as low as —0.6. 

It is necessary to make a theoretical prediction in SM as precision as we can in order to give a firm ground for 
finding new physics. For the purpose, we calculate the twist-3 and weak annihilation contributions in SM using the 
method in Ref. '6II by which there is no any phenomenological parameter introduced. The numerical results show 
that the annihilation contributions to Br are negligible, the twist-3 contributions to Br are also very small, smaller 
than one percent, and both the annihilation and twist-3 contributions to the time-dependent CP asymmetry are 
negligible. The conclusion remains in SUSY models and consequently we neglect the annihilation contributions in 
numerical calculations. 

In conclusion, we have shown that the recent experimental measurements on the time-dependent CP asymmetry in 
B 4>Ks, which can not be explained in SM, can be explained in the SUSY SO(IO) grand unification theories where 
there are flavor non-diagonal right-handed down squark mass matrix elements of second and third generations whose 
size satisfles all relevant constraints from known experiments (r fj.j, B B^ —f fj,'^ fj,~ , B — > Xgfi^ fi~ , B — s- 

X^g, AMs, etc.). Therefore, if the present experimental results remain in the future, it will signal the significant 
breakdown of the standard model and that the SUSY SO(IO) GUT is a possible candidate of new physics. 
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Appendix A Loop functions 

In this Appendix, we present the one-loop function of Wilson coefficients in this work. 

,9/2,12,4,34,6,bo(a;) 



-F2, 12, 4, 34, 6, 60(2;) — X 
^2,12,4,34,6,bo(''') 



dx 

x"^ 12, 4,34,6, 60(2:) 



2 dx"^ 
with 

/i2(x) = 9fi{x) + /2(x), h^ix) = 9/3(0;) + Uix)/2 

where /i,2,3,4,6,bo are given in Ref. p^ . 



